A link between the spin fluctuation and Fermi surface in high Tc cuprates 
— A consistent description within the single-band Hubbard model 



o\ ■ 
a^ ■ 
a^ : 
^ . 

in ■ 

<\ 

\o ■ 



Kazuhiko Kuroki, Ryotaro Arita, and Hideo Aoki 
Department of Physics, University of Tokyo, Hongo, Tokyo 113-0033, Japan 

(February 1, 2008) 

A link between the spin fluctuation and the "fermiology" is explored for the single-band Hubbard 
model within the fluctuation exchange (FLEX) approximation . We show that the experimentally 
observed peak position of the spin structure in the high Tc cuprates can be understood from the 
model that reproduces the experimentally observed Fermi surface. In particular, both the variation 
of the incommensurability of the peak in the spin structure and the evolution of the Fermi surface 
with hole doping in La2-2:Sra;Cu04 may be understood with a second nearest neighbor hopping 
decreasing with hole doping. 

PACS numbers: 71.10.Fd, 74.25.Ha, 74.72.-h 
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In the physics of high Tc cuprates, it is of great im- 
portance to investigate how the states dominated by the 
electron-electron interaction are affected by the band 
structure and the band fiUing. This is especially so for 
the spin fluctuation, which has been recognized as one 
of the central normal state properties. So it is an im- 
portant challenge to look into the link between the spin 
fluctuation and the "fermiology" . 

The doping dependence of the spin fluctuation in 
La2-a;Sra;Cu04 (LSCO) has been well established from 
neutron scattering experiments, where incommensurate 
spin structures appear when holes are doped. It has 
also been shown that such an incommensurate spin struc- 
ture can be understood qualitatively within the single- 
band Hubbard or t-J models. [^-D However, the doping 
dependence of the incommensurability 5 (the deviation 
of the peak in the spin structure from q — (tt, tt)) has 
not been, to our knowledge, understood quantitatively 
within single-band models. Namely, large 5 appears with 
a small amount of hole doping (e.g., 5 ~ 0.2 already at 
X ~ 0.1) in the underdoped regime, while 5 saturates to 
0.25 ^ 0.30 in the overdoped regime. |^ 

More recently, the incommensurate spin structure is 
observed in the (underdoped) YBa2Cu306.6 (YBCOg.e) 
and (nearly optimally doped) Bi2Sr2CaCu208 (Bi2212) 
samples as well, JTo|-p^ which suggest that the phe- 
nomenon is not an accident for LSCO. Several authors 
have theoretically investigated the Hubbard and t-J 
models in the context of these materials. By con- 
trast, it is well known that the antiferromagnetic order 
is robust against electron doping in Nd2-a;Ce2:Cu04_y 
(NCCO). The difference is endorsed by a recent observa- 
tion of a commensurate spin fluctuation in a supercon- 
ducting NCCO sample [x = 0.15). (l| 

On the other hand, the shape of the Fermi surface 
has been revealed from the angle resolved photoemission 
(ARPES) studies for various high Tc cuprates including 
YBCO |5|, Bi2212 (l6|jl|, LSCO 019), and NCCO 
p0|. In particular, a recent ARPES experiment [|18|Jl9| 



has shown that the nesting of the Fermi surface in the 
underdoped LSCO is not so good as has been expected 
previously (2l]j2|] from band calculations. 

So the theoretical challenge is the following. We can 
simulate these shapes of the Fermi surface by taking ap- 
propriate values for the hopping parameters extended to 
distant neighbors in the tight-binding model, and the 
question is how the electron-electron interaction would 
dictate the spin structure. This is exactly the purpose of 
the present study. We employ here the single-band Hub- 
bard model, one of the simplest models for repulsively 
interacting electrons. We shall show that the experimen- 
tally observed Fermi surface and the peak position of 
the spin structure can be understood consistently within 
the single band Hubbard model by taking appropriate 
values for the hopping parameters. The spin suscepti- 
bility and the Fermi surface of the Hubbard model are 
calculated with the fluctuation exchange (FLEX) approx- 
imation, which is appropriate for treating large spin 
fluctuations. 

We consider the single band Hubbard model. 
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H.c) + [/^ 
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where we take tij — t for nearest neighbors, tij = t' for 
second nearest neighbors, and tij — t" for third near- 
est neighbors. Using FLEX approximation, we calculate 
the RPA-type spin susceptibility w) as well as the 
Fermi surface. In the actual calculation, we take 64 x 64 
fc-point meshes and the Matsubara frequencies w„ from 
-(2A^c - l)7rT to (27Vc - 1)ttT, where = 1024. The 
dynamical spin susceptibility for cj ^ is obtained by an 
analytical continuation of x(q, itOn) using Pade approx- 
imation. |2^ For all the case studied, we take U — At 
and T — 0.05t. We denote the band filling as n, the hole 
doping level as ph ~ 1 ^ n, and the electron doping level 
as Pe = ^ — 1- 
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FIG. 1. The Fermi surface for n = 0.85 (a) or n = 1.15 (b) 
with U = At, t' = -0.40t, t" = Omt, T = 0.05t is plotted by 
gray circles. The experimentally observed Fermi surfaces of 
YBCO (Q (a)) lH, Bi2212 (□ in (a)) 0, and NCCO (Q 
in (b)) pd] are superposed for comparison. 



First let us look at the case of YBCO, Bi2212, and 
NCCO. We find that the Fermi surface obtained by 
ARPES experiments for these materials can be repro- 
duced by fitting t' = -0.40, t" = +0.05. These values 
are close to the ones adopted in some of the previous 
studies. p|,|6|,p^ 22 2^ In Fig.|l], we superpose the experi- 
mentally obtained Fermi surface to that of the Hubbard 
model obtained with FLEX. The fit for YBCO is not so 
good as the others because the Fermi surface splits due 
to the interlayer coupling. 

Fixing (t',t") at (-0.40,0.05), we have calculated the 
spin susceptibility xIQj'^) for various band fillings. In 
Fig.^, we show typical results for the imaginary part 
(= x") of xili^) with uj ~ 0.05t for hole-doped (a) or 
electron-doped (b) cases. We can see that incommen- 
surate peaks appear at q = (tt, 7r(l ± 6)), (7r(l ± S),tt) 
for hole doping, while a commensurate peak appears at 
q = (tt, tt) for the electron doping. 

The asymmetry between the hole and electron dop- 
ing becomes even clearer if we plot the incommensura- 
bility (5 as a function of the band filling in FigJ|. 
One can immediately notice that a rather large incom- 
mensurability appears with a small amount of hole dop- 
ing {S ~ 0.25 at ph ~ 0.1), while the commensurate 
spin structure is robust against the electron doping up to 
Pe — 0.25. The robustness of the commensurate spin cor- 
relation in the electron doped regime is consistent with 
previous numerical studies for the Hubbard model |^] 
and the t-J model. |^,^ We have further found that 
the obtained onset of the incommensurate spin structure 
around n ~ 1.3(pe ~ 0.3) in the electron doped regime 
coincides with the situation at which the Fermi surface 
becomes too bloated to intersect the magnetic Brillouin 



zone boundary (the line | fc^^ | -I- | fcy | = tt), as shown in the 
inset of Fig.|. |§ 

The result agrees quantitatively with experimental re- 
sults for YBCOe.e, Bi2212, and NCCO, plotted as open 
symbols in Fig.^. Namely, YBCOg.e, in which the hole 
doping level in the Cu02 planes is around ph 0.1, js^] 
is shown to have 6 ~ 0.22. ||ll|,|l2| For Bi2212, a neutron 
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FIG. 2. x"i<i,^ = 0.05i) plotted against q forn = 0.90 (a) 
or n = 1.15 (b) with U = 4t, t' = -0.40t, t" = 0.05t, and 
T = 0.05t. 
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FIG. 3. Incommensurability S plotted as a function of band 
filling n for [/ = At, t' = -OAOt, t" = 0.05i, T = 0.05t 
(solid squares). We plot here 5 defined for xi^j^ ~ 0), which 
coincides with S defined for x"(cii'^ 0) within 0.02 for 
the cases studied. Also plotted are the experimental results 
for YBC06.6 (O), Bi2212(D), and Ndi.s5Ceo.i5Cu04-a (A). 
The inset shows the Fermi surface for n = 1.30. 



scattering result has been obtained for a nearly optimally 
doped sample, in which an incommensurate spin struc- 
ture has been detected. Although the actual spatial pat- 
tern of the spin structure has not been determined pre- 
cisely, it has been concluded in Ref. ||T^ that 6 would 
be around 0.32 if the spatial symmetry of the spin struc- 
ture is the same as that in YBCO. The optimal hole 
concentration for Bi2212 is around ph ^ 0.18 according 
to Ref. pi, which is adopted in the plot. For NCCO, 



the commensurate antiferromagnetic order persists up to 
about X ~ 0.13. Quite recently, a commensurate spin 
fluctuation has been found in the superconducting sam- 
ples Ndi.85Ceo.i5Cu04_y. ||lj] Here we plot this result by 
nominally taking n = 1.15, although the actual doping 
level may depend on y. 

A few comments are in order at this point. Some 
previous studies have adopted a smaller — i'(= 0.17i) 
and a larger t"(= 0.20i)for YBCO. Although this 
parameter set approximately reproduces the experimen- 
tally obtained Fermi surface as well, |^,^ it only gives 
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a broad commensurate spin structure in the underdoped 
regime {ph < 0.1), which crosses over to an incommen- 
surate structure with peaks in the diagonal direction 
((7r(l ± (5), 7r(l ± S))) rather than for (tt, 7r(l ± S)) in the 
overdoped regime {ph > 0.2). 

In fact, we can reduce the value of —t' without much 
changing the shape of the Fermi surface by taking larger 
values of t" , but with a slight reduction of —t' to 
t' = — O.SOi (with t" = O.lOi), the incommensurate 
spin structure at n = 0.90, although still present, is al- 
ready substantially degraded as compared with that for 
t' ~ —OAOt. Thus, among the sets of parameters that 
reproduce the Fermi surface of YBCO, the ones with 
—t' ^ t" seem to be appropriate from the viewpoint 
of the spin structure. Q 

In Ref. 0, a relatively smaU -t'{= 0.25t with t" = 0) 
has been concluded on the ground that a quantum Monte 
Carlo calculation shows an appearance of the peak in the 
spin structure in the diagonal direction for —t' > O.M. 
Wc have also found that the diagonal peaks can appear 
for the present choice of -t' = OAOt (with t" = +0.05t), 
but that occurs only for band fillings n = 0.70 (for which 
the calculation in Ref. was performed) or less, so that 
there is no inconsistency. However, a smaller —t' is not 
acceptable in our view, since, not only the large curva- 
ture of the Fermi surface requires a large —t\ but also 
the appearance of large S with a small amount of hole 
doping requires a large —t' as we shall see in Figjs] be- 
low. Since S is as large as ~ 0.22 for YBCOg.e (with ph 
being as small as 0.1), a large value of ~t' such as the 
one adopted here is considered to be appropriate for this 
material. 

Let us now turn to La2-a;Sra;Cu04, which provides 
an excellent test-bed for studying the link between the 
Fermi surface and spin structure, since (i) the material 
is unique in that the Fermi surface changes its topology 
(i.e., connectivity) between the underdoped {x = 0.1) 
and overdoped {x — 0.3) regimes, as has been observed 
by a recent ARPES experiment, [|l8yi9[ (ii) the over- 
all doping dependence of the spin structure is well es- 
tablished from neutron scattering experiments. First 
let us look into the Fermi surface. The ARPES result 
suggests that the nesting of the Fermi surface for 
x = 0.1 is not so good as has been expected previously 
from band calculations, so one needs to take relatively 
large values for — i' to reproduce that. We can then take 
t' = -0.35^, t" = +Omt, and n = 0.90 to obtain a 
good fit to the experimentally obtained Fermi surface as 
shown in Fig.^(a). On the other hand, the ARPES data 
a.t X — 0.3 cannot be reproduced, as seen in Fig.^(b), 
for fixed {t',t"){= (-0.35<, +0.050) with n = 0.70. In 
fact, we have to assume —t' to be as small as O.lOt for 
a fixed t"{^ OMt), or t' = -0A5t when t" is varied to 
zero, to reproduce the ARPES data. These values are 
actually close to the ones adopted for LSCO in many of 
the previous studies. [|-U|j2l|,|2| We show in FigJ|(d) 
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FIG. 4. The Fermi surface for U = 4t, T ^ 0.05t, 
t" = 0.05t, and (t',n) = (-0.35t, 0.90) (a), (-0.354,0.70) (b), 
(-0. lot, 0.90) (c), (-0. lot, 0.70) (d). The experimentally ob- 
served Fermi surfaces of LSCO for a: = 0.1 (a,c) and x = 0.3 
(b,d) [[L9| are superposed for comparison. 



how the Fermi surface is well reproduced for t' ~ 
-OAOt, t" = Omt with n = 0.70. A change in t' is 
really required, since if we fix {t',t") to these values to 
go back to n — 0.90 in Fig.^(c), the Fermi surface is seen 
to deviate from the ARPES data for a; = 0.1. 

One might suspect that the reason why we have to 
take different values for the hopping parameters in the 
underdoped and overdoped regimes might be because the 
FLEX approximation did not take into account properly 
the electron correlation effect, which should become more 
significant in the underdoped regime. This, however, can- 
not be the case, since it is known that the electron cor- 
relation effect tends to enhance the nesting in Fermi sur- 
faces as the system approaches half- filling, while the ex- 
perimental data suggest that the nesting becomes worse 
for smaller x. Thus we can assume that the decrease of 
—t'/t with hole doping is not an artefact. 

Keeping the above analysis for the Fermi surface in 
mind, we have calculated x(1j'^) for three values of t' , 
namely, t'/t ^ -0.35, -0.25, and -0.10, fixing t" 
omt. A typical result for t' = -0.25t, t" = +0.05t, 
and n = 0.80 is given in the inset of Fig.|[ Again, 
incommensurate peaks are found at k = (tt, 7r(l ± 5)), 
(7r(l ± (5), tt). The band fiUing dependence of the incom- 
mensurability is plotted in Figj^. Here, a result similar to 
that for (t', t") = (-O.lOi, 0.05t) can also be obtained for 
{t',t") = (—0.15^,0.0), i.e., the other choice of param- 
eter set that reproduces the ARPES data for x = 0.3. 
The doping dependence of 6 for small values of —t' is 
similar to those obtained in previous theoretical studies 
for LSCO. |^,|| As already stated previously, |||,|^ S is 
essentially zero in the underdoped regime {n > 0.90) for 
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FIG. 5. A plot similar to Fig.| for U = 4t, T = 0.05t, 
t" = O.OSt, and t' = -0.35t (solid squares), -0.25t (solid cir- 
cles), —O.lOt (solid triangles). The hatched area represents 
the experimental results for LSCO. Q Note S defined here 
is larger by a factor of 2 than that defined in Ref. j^. The 
inset shows a plot similar to Fig.^ for U = 4t, t' = —0.25t, 
t" = 0.05t, T = Omt.Lj = Omt, and n = 0.80. 



small values of — t', which is inconsistent with experi- 
mental observations. It can be seen that the appear- 
ance of large S with a small amount of hole doping , 
which is observed experimentally {S ^ 0.2 for x ^ 0.1), 
1^ can only be obtained for large enough — On the 
other hand, the incommensurability increases monotoni- 
cally with hole doping if we stick to a fixed {t' , t"). Thus 
we conclude that we have to have —t'/t (and possibly 
t" /t) decreasing with hole doping in order to obtain the 
experimentally observed saturation of 5. |^] The doping 
dependence of both the Fermi surface and the spin struc- 
ture can then be understood within a single scenario. 

To conclude, the main message of the present paper is 
that the experimental observations such as the saturation 
of the incommensurability in LSCO and the sharp con- 
trast of the spin structure between electron doped and 
hole doped systems can be understood within a single- 
band model, and that we do not have to go back, as one 
might suspect, to the original three-band model that ex- 
plicitly takes into account Cu 3c? and O 2p orbitals. 

A remaining problem is why, if our picture is correct, 
—t'/t (and possibly t" /t) decreases with hole doping in 
LSCO. As for the discrepancy between the ARPES result 
and the band calculation in the underdoped regime, we 
may make a following conjecture at the present stage. If 
we closely look at the band calculation for LSCO, [ p3|j24| 
there is, in contrast to that of YBCO, |34| a large dis- 
persion between F (k = (0,0)) and Gi (k = (7i',0)) in 
the band lying near the Fermi level. In a tight-binding 
model, this corresponds to a small —t'. A large energy 
difference between F and Gi is known to be due to the 
mixing between the Cu d^2_y2- O Px,y orbital and the 

If, for some reason. 



regime, so that the main character of the band near the 
Fermi level is Cu dx2_y2 - O Px,y, then the energy dif- 
ference between F-Gi would be smaller, amounting to 
a larger —t'. If this is correct, the reason why the hy- 
bridization is weak in the underdoped regime, or why it 
is recovered in the overdoped regime becomes a future 
problem. Experimental observation in Ref. | |35[ | showing 
that the Cu d^^^-r^ - O Pz character grows with hole 
doping may have some relevance to this point. 

The relation between the present result and the so- 
called '1/8-problem' or 'charge stripes' is another 
interesting future problem. 
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